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We have performed coUinear and non coUinear calculations on neutral Bi4,Mn, and coUinear ones 
on ionized Bi^Mn with charges +1 and -1 to find out why theoretical calculations will not predict 
the magnetic state found in the experiment. We have used the density functional theory to find a 
fit between the theoretical prediction of the magnetic moment with the experimental value. Our 
calculations have consisted in a structure search of local energy minima, and then a search of the 
magnetic lowest energy state for each resulting isomer. The geometry optimization found 3 local 
minima whose fundamental state is the doublet spin state, which could not be found in previous 
theoretical works, but they are higher in energy than the lowest-lying isomer by « 1.75 eV. This 
magnetic state could help understand the experiment. Calculations of non-coUinear magnetic states 
for the BiiMn do not lower the total magnetic moment. We conclude arguing how the 3 isomers 
with doublet state could actually be the ones measured in the experiment. 



I. INTRODUCTION 

Nanostructured materials have become the main 
source in current technology to develop devices with 
novel properties, as low dimensionality causes matter to 
behave differently than the known properties of bulk. 
Among them, the zero-dimensional materials, i.e. atomic 
clusters and nanoparticles, have attracted a great deal of 
attention, specially their magnetism. Pure small clusters 
can exhibit far larger magnetic moment per atom than 
the isolated atom or their crystals [l|-l3|, and this effect 
can be enhanced when impurities are added [4'|. Within 
this trend, the study of transition metal binary clusters 
was triggered by the exceptionally high magnetism found 
in CoRh nanoparticles by Zitoun et al.;f5| since then a 
number of magnetically enhanced nanoalloys of ferromag- 
netic and non-magnetic transition metals have been stud- 
ied theoretically [6|-l9| and experimentally [10|-ll2|. Yin et 
al. |13| found an enhanced magnetism in BiMn clusters for 
Bi-to-Mn ratios close to 2 in their Stern-Gerlach mea- 
surements, which deviated from the trend in other Bi 
nanoalloys like BiCo, where the magnetic moment had a 
small dependence with cluster size[10|. In a later paper, 
Chen et al.|14'] perform extensive density function the- 
ory (DFT) calculations on BinMrim (n=l-6, m=l-12) 
to learn their structure and how their magnetism works. 

While their calculations have been found to be in quite 
a good agreement with the experiment for this cluster 
series, there are a few instances of large discrepancies 
between the theoretical value of the magnetic moment 
and the measured one. The most significative one, as 
highlighted by Chen et al., is the Bi^Mn where DFT 
computation predicts a total magnetic moment, fiT — 
5fiB while the experiment measures 1.6/iB. Extensive 



geometric optimization and an estimation of the orbital 
contribution to the magnetic moment were performed to 
solve this problem, but the cause of this discrepancy was 
not found. 

Addressing the reason for these disparities is key. One 
one hand, the existence of an unadressed question in the 
literature calls for it to be explained. On the other hand, 
these problems might point out defficiencies in our ap- 
proach to ab initio calculations, or in the experimen- 
tal techniques used. The relevance of DFT calculations 
comes from the reliability of their predictions, and we 
need to learn the conditions for which they work, and if 
they do not work for a case, we need to learn the reasons. 

In the present paper we will explore the possible 
sources of the mismatch between theoretical calcula- 
tions and experimental results for the BiiMn nanoalloy, 
through DFT calculations. The considered possibilities 
are: finding isomers of the cluster whose total magnetic 
moment are low enough; allowing non-collinear magnetic 
states for the studied BiiMn structures, as this might 
lower their magnetic moment; finding out the orbital con- 
tribution to the magnetic moment of each structure. We 
will organize our communication as follows. In Sec. II 
we explain the computational methods we have used to 
perform the different types of calculations we performed, 
also reasoning why each method was used seeking which 
result. In Sec. Ill we show our results and discuss them, 
this section being divided in two subsections: subsec. 
1 contains the results of the calculations with imposed 
coUinear magnetic moments, and subsec. 2, the results 
of the calculations where the coUinear is lifted, allowing 
magnetic moments to arrange freely. Sec. IV contains 
the conclusions derived from our findings. 



II. COMPUTATIONAL DETAILS 



III. RESULTS 



Our first step was to perforin ab initio calculations 
within the framework of the DFT over a large number 
of different structures of Bi/^Mn, both taken from the 
literature and created by the authors. We have solved 
each system using linear combinations of Gaussian-type 
orbitals within the Kohn-Sham (KS) density functional 
methodology (LCGTO-KSDFM), with the deMon 2003 
code 15^. The calculation of the exchange-correlation 
(XC) energy term was carried out using the generalized 
gradient approximation (GGA), in the Perdew-Burke- 
Ernzerhof ansatz|16|. It has been shown that the GGA 
is the approximation that gives the best results for ge- 
ometry optimizations of small cluster s 1171 KLSf . The or- 
bital basis set used in our calculations[19'| for Bi and Mn 
uses the effective core potential from Stuttgart-Dresden: 
RECP — SD[20|. It considers a relativistic approxima- 
tion for the outer 15 electrons of the Mn, and a quasi- 
relativistic one for the 5 outer electrons of the Bi. The 
electron density is expanded in auxiliary basis sets in or- 
der to avoid the calculation of the N'^ scaling Coulomb 
repulsion energy, where N is the number of the basis 
functions. The auxiliary basis sets are GEN- A3* 
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the Bi and A2-DZVP [22] for Mn. Each structure was 



geometrically optimized through a Born-Oppenheimer 
molecular dynamics simulation. After this first approach, 
the geometry of each cluster was optimized further with a 
Broyden-Fletcher-Goldfarb-Shanno algorithm [2 3]. Con- 
secutively, we carried out a search in the collinear mag- 
netic state that yields the energy minimum for zero ex- 
ternal magnetic field and T=0 for each geometry. After 
this optimization, an initial structure to which we impose 
different magnetic configurations often relax to different 
final structures, while different initial structures might 
converge to the same final structure. Henceforth we ob- 
tain a large number of relaxed isomers; we focus on the 
21 geometries with lower energy. At this point we per- 
form single-point calculations for each isomer, searching 
for their fundamental magnetic state. Next we proceed 
to perform non-collinear single-point calculations to de- 
termine the lowest-energy magnetic state for each one of 
the 21 isomers. These calculations are performed using 
a different DFT implementation, the Octopus code[2J]. 
The Octopus calculations have considered wavefunctions 
to be complex spinors and they include the spin-orbit 
coupling relativistic correction. In the same manner as 
the collinear calculations, XC energy is approximated 
through a GGA-PBE functional. The inner electrons 
of the atoms are represented with a core model, u sing 
the Hartwigsen-Goedecker-Hutter pseudopotentials.[25J 
As for the grid upon which the KS equations are solved, 
we have chosen a radius of 4.4 and a spacing of 0.1. The 
mixing scheme is a broyden one, with a mixing factor of 
0.02. 



1. Collinear Calculation Results 

In Fig. [1] we represent the 21 lowest-energy isomers 
that we have found in our calculations and in Fig. [2] we 
plot the energies of each of their magnetic states from 1 
to 9/iB. Our lowest energy structure, an edge-capped dis- 
torted tetrahedron, and its ^itotai — ^I^^b agree with the 
previous theoretical calculations, which in turn disagree 
with the experimental data[l3|. Also, we have found the 
pyramidal structure to be 0.76 eV higher than the lowest 
lying isomer, and the W-like structure 1.06 eV higher, 
which is in very good agreement with Chen's work. 14 
Furthermore, isomers with low coordination for the Mn 
tend to have lower energies as well. 

The structures with fj-totai = 5fj,B dominate the lowest 
energy isomers until the 11*'' one, and from that point 
on, structures with fitotai ~ SfiB dominate instead until 
the the 18*'' one. Our 3 highest energy isomers have 
^J'total = 1/^B- Studying separately the different isomers 
for each of these 3 magnetic states, we have found that 
the energy of clusters set to Utotai — '^P^b and 5/iB grows 
in average with the ground state energy. On the other 
hand, if the clusters are kept at fitotai = l/^s, we find no 
such tendency. 

Analysing each atom of the nanoalloy separately, most 
of the magnetic moment comes from the Mn atom, while 
the Bi magnetic moments tipically remain close to zero. 
As we can see in Fig. |3K), electrons are shared from the 
Bi to the Mn, in particular from 6p orbitals to Ad. This 
transfer does not generally enhance magnetism in either 
element, as the Mn total magnetic moment oscillates be- 
tween 4 and S/is for all clusters, without correlation to 
charge transfers (see Fig. [3b)). In the 3 higher energy 
isomers, however, one of the Bi atoms stands isolated, 
bound only to the central Mn. This Bi, in addition to 
transferring 6p electrons to the Mn, also receives charge 
into its 7p orbital, greatly increasing its magnetic mo- 
ment which couples antiparallel with the Mn. This makes 
the total magnetic moment of the isomer to decrease to 
IfiB- The lower total magnetic moment state is hence 
originated in the enhancement of the magnetism of the 
Bi atoms, whose magnetic moments are either very close 
to zero or antiparallely aligned to those of the Mn atoms. 

We have not found the lowest energy isomer to have a 
total magnetic moment close to the experimental value. 
We have found the structures whose ground state is close 
to the experimental total magnetic moment to be ener- 
getically close enough to the lowest energy one to have a 
large enough population in the synthesis to justify the ex- 
perimental results. So we proceed to test other exchange- 
correlation functionals to find out if this could be the 
reason for the discrepancy. In Fig. 2] we show the ground 
states calculated for the lowest-lying isomer for different 



exchange-correlation (XC) functionals, and we compare 
them with our results for PBE functional. We performed 
this calculation for the values of ht = 1,3, 5^b. Lo- 
cal functionals exhibit higher energies, but all functionals 
provide the same energetic ordering of the three magnetic 
states: fir — ^IJ-b as the lowest, 3/LtB second lower and 
1/iB third lower. When splitting the energy in XC and 
classical terms, all functionals keep the same ordering: 
XC term greatly favors the 3/zb and 5/is states, while 
the bias of the classical term lowering the 1/is state can- 
not balance it. 

The final scenario we can examine with coUinear cal- 
culations is the possibility that the experimental BiMn 
clusters are ionized in the moment of their production, 
or that their structure changes during the ionization in 
Stern-Gerlach. This is unlikely because of the low tem- 
perature of the process. [26[ The results still do not yield 
any significative reduction of the total magnetic moment, 
as we can see in Fig. [5l The lowest lying isomers are small 
deformations of the one for the neutral cluster, and the 
lowest magnetic configuration is for either ion 4/xb. The 
lowest magnetic state is found for the anion third lowest 
lying isomer, which has fiT = ^fJ-B- Even this result is 
larger than the experimental value, so it cannot explain 
the discrepancy. 

As a summary, we have found that there are 3 struc- 
tures of the Bi^Mn cluster whose fiT = IfJ-B- The ex- 
perimental value is 1.6/xb, which implies that it is not 
isomerically pure: there has to be a population mixture 
of clusters with different structures and different fiT , in- 
cluding IfiB- But these 3 structures are « l.TSeF higher 
than the lowest-lying one, and there are 18 structures 
with lower energy than them. Consequently, it is un- 
likely that the experimental device can produce these 
isomers in a high enough ratio as to lower the magnetic 
moment, explaining the experimental value. Except that 
some other reason, unaccounted for in our calculations, 
happens to favour their synthesis. 



2. Non-Collinear Calculation Results 

The next step is finding out if the discrepancy can be 
solved with non-collinear calculations, and we perform 
such calculations for the 21 lowest-energy structures we 
found in the collinear calculations. In table U we show 
the results, with the energy and the components of the 
total magnetic moment for each geometry. As we can see 
in Fig. |B^), the energy dependence of the different iso- 
mers with non-collinear magnetic configuration roughly 
follows that of the isomers with collinear magnetism. The 
total magnetic moment (see Fig. [6j))) ranges from 3.7 to 
6.4^B, with the lowest lying isomer having /j,t = 5.0/iB, 
values which exceed the experimental value and agree 
with the collinear calculations for the lowest energy clus- 
ters. In the last column of the Table U we compare the 



modulus of the total magnetic moment components: M^ 
and My against M^ . This ratio is very small, specially for 
the lower energy isomers, hence we can conclude that the 
collinear approximation is a good one, which validates 
the results of our collinear calculations in the previous 
subsection. 

Further analysis of the magnetic components of our 
isomers yield the orbital magnetic moment: In Fig. [BJ:) 
we compare the modulus of the orbital magnetic moment 
with the value for BiMn in bulk, 0.17/j,b, 27[ seeing that 
the value of the orbital magnetic moment has the same 
order of magnitude, being small in all cases. In the col- 
umn 5 of TableUwe compare the orbital component with 
the total magnetic moment of each isomer, and we find 
that in all cases, this contribution is smaller than 10%, 
except for the structure 16, where the ratio between the 
moduli is 14%. From our results, we conclude that the 
orbital magnetic moment for these clusters is too small 
to be a significative factor in the total magnetic moment, 
so it cannot be the source of the discrepancy between 
the experimental results and the theoretical calculations. 
Furthermore, our calculated orbital magnetic moments 
for each isomer are well within the upper limit estimated 
by Chen et al., 1.35/iB. 



IV. CONCLUSIONS 

We have explored the possible scenarios that could 
lead to the known discrepancy between the experimental 
value of the total magnetic moment of the BiiMn cluster 
and the theoretically calculated one. We have found that 
three isomers among the 21 most stable actually have a 
total magnetic moment below the experimental one, but 
these isomers are too high in energy respect to the lowest 
lying one. Furthermore, there are 18 isomers with lower 
energy, so these three cannot make up for a fraction of 
the population of randomly created Bi^Mn clusters sig- 
nificative enough to reduce the average total magnetic 
moment down to the experimental value. We have not 
found either that using different XC functional approxi- 
mations makes the ground state of the lowest lying isomer 
have a lower magnetic moment. The analysis of the pos- 
itive and negative singly-charged ions has also yielded no 
lowest energy structure with a magnetic moment closer 
to the experimental one. In fact, none of the ions of the 
21 isomers we have considered has a ground state with 
a magnetic moment lower than the experiment. This 
should be expected, though, as the only lower magnetic 
state available to ions of Bi^Mn is QjjLb- The possible 
presence of ions would not help us explain the experi- 
ment unless we had found ground states at /^t = O/is. 
From our non-collinear calculation, we rule out other two 
possible sources of the disagreement: the magnetic con- 
figuration is collinear to a high degree, even when allow- 
ing the individual atomic magnetic moment to arrange 
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TABLE I: Spread of the results obtained in our non-collinear 
calculations. First column labels each structure according to 
their energy ordering in the collinear calculations. Second 
column units are electronvolts. Column 3 is the modulus of 
the total magnetic moment and units are /is. Column 4 is the 
modulus of the orbital moment and units are ^b ■ Column 5 
is the dimensionless ratio of the orbital moment and the total 
magnetic moment. Column 6 is the dimensionless ratio of the 
XY plane component of the total magnetic moment and the 
Z axis component. 



structure 19, 20, 21 are very similar the lowest lying iso- 
mer for Bi^Mn clusters (a tetrahedron), simply adding 
an extra Bi atom to the Mn end of the tetrahedron. If 
Bi^Mn clusters form much more quickly in the experi- 
mental device than Bi^Mn ones, it is not unreasonable 
to think that the later will form from the former. Bi^Mn 
also has a large electrostatic dipole we have calculated to 
be 2.4 D in the direction that connects the Mn atom with 
the center of the triangle formed by the three Bi atoms. 
This dipolar moment could help a fourth Bi atom to cou- 
ple to the Mn instead of breaking up the tetrahedron to 
form the calculated Bi^Mn lowest lying structure. If the 
process of measurement of mass and magnetic moment of 
the clusters is fast enough, they might not have enough 
time to relax into the said lowest lying isomer, hence re- 
sulting in some of our structures 19, 20, 21 making up a 
large fraction of the measured Bi^Mn clusters. 

To verify this, one possible path would be to perform 
an analysis of the optical properties of the experimen- 
tally obtained of the BiiMn clusters to identify their 
geometry, and see if they match with the isomers theo- 
retically obtained as the lowest energy ones, or instead 
they match with any of the 19, 20, 21 isomers that actu- 
ally show a low total magnetic moment. If the later case 
were true, then we would have to explain why the ex- 
perimental setup produces structures that theoretically 
at T=0K are known not to be the fundamental one. On 
the other case, it would be necessary to know why the 
theory cannot predict the correct magnetic moment for 
a cluster with known composition and structure. 
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FIG. 1: (color online) Geometry of the 21 lowest energy struc- 
tures for the nanoalloy BiiMn, from 1) lowest to 21) higher. 



1.2 



1.0- 



0.8- 



O 0.6- 
CD 

I 

^0.4^ 



0.2- 



0.0- 




"T" 
9 



15 



18 



21 



12 

structure 



FIG. 2: (color online) Energy comparison of the 21 structures 
with lower energy for the magnetic configurations utotai ~ 
1, 3, 5, 7, 9/is- The numbers in the x-axis correspond to the 
structures shown in Fig. [T] 
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FIG. 3; (color online) Atomic charges (a) and magnetic mo- 
ments (b) of each atom of the ground state of each one of the 
21 lowest energy isomers of Bi/^Mn. 
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FIG. 4: (color online) Comparison of the 3 lowest magnetic 
states of the lowest lying isomer of BiiMn calculated using 
different XC functionals. 
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FIG. 5: (color online)Energy comparison of the 21 struc- 
tures with lower energy for the magnetic configurations /it ~ 
0, 2, 4, 6, 8/iB of the negatively (a) and positively (b) charged 
ions of BiiMn. 
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FIG. 6: (color onlinG)Encrgy (a), magnetic moments (b) and 
orbital moments (c) for the 21 lowest energy isomer oi BiiMn 
found in the collinear calculations. In c), the horizontal line 
marks the value of the orbital magnetic moment in bulk BiMn. 



